B^^-CONF-07/012 
SLAC-PUB-12734 
arXiv:0708. 1630 [hep-ex] 
August 2007 



Time-Dependent Analysis of the Decay p^p^ 

The BABAR Collaboration 
August 12, 2007 

Abstract 

We study the decay — > p^p^ in a sample of about 427 million T{AS) — > BB decays collected 
with the BABAR detector at the PEP-II asymmetric-energy e~^e~ collider at SLAC. We find the 
branching fraction B = (0.84 ib 0.29 it 0.17) x 10~^ and longitudinal polarization fraction of = 
0.70 lb 0.14 lb 0.05, where the first uncertainty is statistical, and the second is systematic. The 
evidence for the B^ — > p^p^ signal has 3.6a significance. We investigate the proper-time dependence 
of the longitudinal component in the decay and measure the CP-violating coefficients = (0.5 ib 
0.9 lb 0.2) and C£° = (0.4 ib 0.9 ib 0.2), corresponding to the sine and cosine terms in the time 
evolution of asymmetry. We study the implication of these results for penguin contributions in 
B ^ pp decays and for the CKM unitarity angle a. 
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1 INTRODUCTION 



Measurements of CP-violating asymmetries in the B^B^ system test the flavor structure of the 
standard model by over-constraining the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix [1]. The time-dependent CP asymmetry in the decays of B^ or B^ mesons to a CP eigenstate 
dominated by the tree- level amplitude b — > uud measures sin 2aes , where aes differs from the CKM 
unitarity triangle angle a = arg [—V^^Vfl/V^^V*^^] by a quantity Aa accounting for the contribu- 
tions from loop (penguin) amplitudes. The value of Aa can be extracted from an analysis of the 
branching fractions of the B decays into the full set of isospin-related channels |:2J. 

Branching fractions and time-dependent CP asymmetries in 5 — > vrvr, pn, and pp have already 
provided information on a. Since the tree contribution to the B^ — > p^p^ [3] decay is color- 
suppressed, the decay rate is sensitive to the penguin amplitude. The B^ p^p^ decay has a much 
smaller branching fraction than B^ — > p'^p~ and B~^ — > p'^p^ channels [H 13 El [3 [8] , which leads 
to a more stringent limit on Aa from isospin analysis [21 \7\ [12] than is possible in tttt system. This 
makes the pp system particularly effective for measuring a. 

The error due to the penguin contribution becomes the dominant uncertainty in the measure- 
ment of a using B ^ pp decays. However, in contrast to -B — > vr'^vr'' decays, the four-track final 
state makes a time-dependent analysis of B^ — > p^p^ decays feasible. It allows us to measure the 
CP parameters and Cj^ directly, analogous to and C^~ , resolving ambiguities inherent 
to isospin triangle orientations. The C^^ coefficient is associated with the difference in decay am- 
plitudes for B — > p^p^ and B — > p^p^, while the S^^ coefficient involves interference between the 
B^-B^ mixing and decay amplitudes. 

In B ^ pp decays the final state is a superposition of CP-odd and CP-even states. An isospin- 
triangle relation [2] holds for each of the three helicity amplitudes, which can be separated through 
an angular analysis. The helicity angles 6i and 02 are defined as the angles between the direction 
of vr^ and the direction of the B in the rest system of each of the candidates. The resulting 
angular distribution is given by 

d^T/{T dcos 01 dcos ^2) = ^ | ^ (1 - fi) sin^ Oi sin^ 02 + fi cos^ 0i cos^ 02 1 , (1) 

where Jl = \Aq\'^ / {T,\Ax\'^) is the longitudinal polarization fraction and ^;^=_i^o,+i are the helic- 
ity amplitudes. The fraction of longitudinal polarization is a priori unknown. Polarization was 
expected to be predominantly longitudinal [9]. However, significant departure from this expec- 
tation was found in penguin-dominated i?-decay modes [TO] and polarization and CP-asymmetry 
measurements in the B^ p^p^ decay may help in resolving this puzzle [llj . 

In this paper, we update our previous measurement of the branching fraction and longitudinal 
polarization fraction in B^ —>■ p^p^ decays |1], and present the first study of the time-dependent 
CP asymmetry Acp in this mode. We determine the coefficients C^" and S^^ of the asymmetry 
for the longitudinal component, which is given by 

Acp{At) = -Cf cos AmAt + 5f sin AmAt (2) 

These coefficients, together with improved measurements of the branching fraction and longitudinal 
polarization, allow a complete isospin analysis and improved constraints on the penguin contribution 
to B ^ pp decays. Changes with respect to our previous analysis [3] include increased datasample, 
improved track-selection techniques, and inclusion of the i?-decay time information. 
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2 DETECTOR AND DATASET 



These results are based on data collected with the BABAR detector [13) at the PEP-II asymmetric- 
energy e~^e~ collider p^. A sample of 427 it 5 million BB pairs was recorded at the T{4S) 
resonance with the center-of-mass (cm.) energy ^/s = 10.58 GeV. Charged-particle momenta and 
trajectories are measured in a tracking system consisting of a five-layer double-sided silicon vertex 
tracker and a 40-layer drift chamber, both within a 1.5-T solenoidal magnetic field. Charged- 
particle identification is provided by measurements of the energy loss in the tracking devices and 
by a ring-imaging Cherenkov detector. 

3 ANALYSIS METHOD 

We select B M1M2 —>■ (7r+7r^)(7r"'"7r~) candidates, with Mi^2 standing for a or /o candidate, 
from neutral combinations of four charged tracks that are consistent with originating from a single 
vertex near the e~^e~ interaction point. We veto tracks that are identified as kaons or electrons. 
The identification of signal B candidates is based on several kinematic variables. The beam-energy- 
substituted mass, TTiES = [(•s/2 + Pi • PB)'^/Ef — p^]-*^/^, where the initial e+e" four-momentum 
(-Bj, Pi) and the B momentum pB are defined in the laboratory frame, is centered near the B mass 
with a resolution of 2.6 MeV/c^ for signal candidates. The difference AE = — -v/i/2 between 
the reconstructed B energy in the cm. frame and its known value has a maximum near zero 

with a resolution of 20 MeV for signal events. Four other kinematic variables describe two possible 
TT+TT" pairs: the invariant masses mi, m2 and the helicity angles 0i, 02- 

The selection requirements for signal candidates are the following: 5.245 < mES < 5.290 GeV/c^, 
\A.E\ < 85 MeV, 550 < rn-1^2 < 1050 MeV/c^, and | cos0i^2| < 0.98. The last requirement removes a 
region corresponding to low-momentum pions with low and uncertain reconstruction efficiency. In 
addition, we veto the copious decays B^ D(*)-7r+ ^ {h+ TT TT )7r^, where refers to a pion or 
kaon, by requiring the invariant mass of the three-particle combination to differ from the Z)-meson 
mass by more than 13.2 MeV/c^, or 40 MeV/c'^ if one of the particles is consistent with a kaon 
hypothesis. 

We reject the dominant e+e~ qq (q = u, d, s, c) (continuum) background by requiring 
I cos I < 0.8, where 6t is the angle between the S-candidate thrust axis and that of the re- 
maining tracks and neutral clusters in the event, calculated in the cm. frame. We further suppress 
continuum background using a neural network-based discriminant £, which combines eight topolog- 
ical variables calculated in the cm. frame. In addition to cos 9t, they are the polar angles of the B 
momentum vector and the -B-candidate thrust axis with respect to the beam axis, the value of the 
event thrust, two Legendre moments Lq and L2 of the energy flow around the -B-candidate thrust 
axis [15] computed separately for neutral and charged particles, and the sum of the transverse 
momenta of all particles in the rest of the event, calculated with respect to the B direction. 

We use multivariate S-flavor-tagging algorithms trained to identify primary leptons, kaons, 
soft pions, and high-momentum charged particles from the other B, called -Btag [12] • The effective 
tagging efficiency, which takes into account the efficiency to find the tag and the mistag probability, 
is (31.1 ±0.3)%, as determined on a sample of fully reconstructed open-charm decays. We use both 
tagged and untagged events in our sample. Additional background discrimination power arises 
from the difference between the tagging efficiencies for signal and background in seven tagging 
categories (ctag = 1..7). We determine the proper time difference At between the signal B and 
-Btag from the spatial separation between their decay vertices. The -Stag vertex is reconstructed 
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from the remaining charged tracks in the event and its uncertainty dominates the At resolution 
aAt- The average proper time resolution is (o"At) ~ 0.7 ps. Only events that satisfy |At| < 15 ps 
and fJAt < 2.5 ps are retained. 

After application of all selection criteria, A'cand = 65637 events form the sample for the max- 
imum likelihood fit. On average, each selected event has 1.05 signal candidates, while in Monte 
Carlo [16j samples of longitudinally and transversely polarized p^p^ decays we find 1.15 and 

1.03 candidates, respectively. When more than one candidate is present in the same event, the 
candidate having the best consistency with a single four-pion vertex is selected. Simulation 
shows that 18% of longitudinally and 4% of transversely polarized — > p^p^ events are misre- 
constructed with one or more tracks not originating from the B^ p^p^ decay. These are mostly 
due to combinatorial background from low-momentum tracks from the other B meson in the event. 
Such events still carry the characteristic topology of B^ — > p^p^ events; they are modeled separately 
from the perfectly reconstructed events and included into the probability density functions. 



We use an unbinned extended maximum likelihood fit to extract the B" — > p"p" event yield and 
fraction of longitudinal polarization fi. We also fit for the event yields of B^ — > p^/o and B^ /o/o 
decays, as well as yields in several background categories. The likelihood function is 



where rij is the unconstrained (except if noted otherwise) number of events for each event type 
j (5° p^p^ , 5° pO/o(980), 5° /o (980) /o (980), several background components from 
exclusive and inclusive B decays, and continuum), and Vj{xi) is the probability density function 
(PDF) of the variables Xi = {niEs, A.E,£, mi, m2, cosOi, cos92,Cta.g, At, a/\t}i for the ith event. 

Since the statistical correlations among most variables are found to be small, we take each 
Vj as the product of the PDFs for the separate variables. In a number of special cases, for 
the mass-helicity PDF of continuum backgrounds, or for the mass and helicity PDFs of the B 
backgrounds, we use 2- or 4-dimensional PDFs to properly describe the kinematic correlation 
between the observables. 

We use double-Gaussian functions to parameterize the mss PDFs for fully-reconstructed signal 
events, double-Gaussian functions for AE and relativistic Breit-Wigner functions for the resonance 
masses of p^ and /o(980) [18]. The angular distribution at production for B^ — > p^p^, B^ — > p^fo, 
and B^ — > /q/o modes (expressed as a function of the longitudinal polarization in Eq. ([1]) for 
B^ —>■ p^p^) is multiplied by a detector acceptance function ^ (cos 6*1, cos ^2)) determined from 
Monte Garlo. The distributions of misreconstructed signal events are parameterized with empirical 
shapes in a way similar to that used for B background discussed below. The neural network 
discriminant £ is described by two (continuum) or three (.B-decay events) asymmetric Gaussian 
functions with different parameters for signal and background distributions. 

The PDFs for inclusive B decay modes are generally modeled with empirical analytical distri- 
butions. Several variables have distributions similar or identical to those for signal, such as rriES 
when all four tracks come from the same B, or vr^vr" invariant mass when both tracks come 
from a p^ meson. Also for some of the modes the two vr^vr^ pairs can have different mass and 
helicity distributions, e.g. when only one of the two combinations comes from a genuine or /o 



4 MAXIMUM LIKELIHOOD FIT 




(3) 
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meson, or when one of the two pairs contains a high-momentum pion (as in i? — > oivr). In such 
cases, we use a four-dimensional mass-helicity PDF. 

The proper-time distribution for signal and background B decays 

e-|At|/r c ^ 
/(At, Q) J— X U - QAw + Qn{l - 2uj) + [Q(l - 2-^) + - QAlo)] Acp{At) \ (4) 

is convolved with a resolution function, modeled as a superposition of three Gaussian distributions 
with the means and widths scaled by the per-event error on At. Here Q is the flavor of -Stag, 
w is the average mistag probability, and Aw and parameters describe the difference in mistag 
probability and the tagging efficiency asymmetry between and B mesons. The time distribution 
of continuum background is assumed to have zero lifetime. 

The signal and i3-background PDF parameters are extracted from simulation. The Monte Carlo 
parameters for m-Es, AE, and £ PDFs are adjusted by comparing data and simulation in control 
channels with similar kinematics and topology, such as B^ D'lr^ with D~ — > K~^7r~7r^. The 
continuum background PDF shapes are extracted from off-resonance data or on-resonance sideband 
data, with parameters of the most discriminating PDFs (mES) AE, £) left free in the final fit. The 
discrete i?-flavor tagging PDFs and parameters of the proper time distributions for signal modes 
are obtained in dedicated fits to events with identified exclusive B decays [17] . The tagging PDFs 
for inclusive B backgrounds are determined by Monte Carlo and their systematic uncertainties are 
studied in data. 

Backgrounds from "charmless" 5 — > n transitions, in particular events containing p, /q, or K* 
mesons, have kinematic distributions similar to those of signal events. We study the contributions 
of the dominant decay modes in high-statistics exclusive Monte Carlo samples. We also develop 
two complementary strategies to model these backgrounds in the likelihood fit. 

In the first approach, we single out contributions from the following dominant modes: B^ — > 
af 7r=F, ^ p^K*^, B^ foK*°, B+ p+ , B+ a?7r+, B rj'K, and B+ p^7r+. The 
contribution of the B^ — > a^vr^ decays includes both events where all four tracks are correctly 
associated to the B candidate, and events where at least one of the tracks is picked from the other 
B meson (so-called B^ afir^ "self-crossfeed" events). The event yield of the B^ — > afn^ decays 
is allowed to vary in the fit, while the yields of other six charmless modes listed above are fixed to 
the expected values [19 ^ 120 1 [21]. The events from open charm b ^ c transitions are parameterized 
as a separate background component, with its yield allowed to vary in the data fit. 

In the second strategy, we split the B background into three distinct categories: B^ — > afn^ de- 
cays where all four tracks are correctly associated with the B candidate, an appropriately weighted 
combination of dominant charmless decays, and the rest of the generic b ^ u, b ^ c, 6— > s, and 
b ^ d transitions. For the charmless event category, we combine the following modes: — > a^^vr^ 
self-crossfeed events, 5° ^ p^K*^, B° foK*°, -8° ^ p+p", B° ^ p^vr^, B+ p+p^, 
B^ ai7r~^, B~^ — > af Jq, B — > r]'K, and B'^ p^ir^. Kinematic distributions in these events, 
especially events in which at least one charged particle was not correctly associated to the B can- 
didate, are similar to each other, and also to other, poorly measured charmless decays. This allows 
us to vary the overall event yield for this category of events, after fixing the relative weights of 
each mode to the expected values [TU [20l |2l]. For yet-unmeasured B{B~^ o-t fo), we assume a 
conservative value of (1 it 1) x 10~^; this branching ratio corresponds to the expectation of 10 it 10 
events in the selected data sample. Event yields associated with B^ — > o^tt^ component and the 
generic B decays are allowed to vary in the maximum likelihood fit. 

We find that both strategies for describing B decay backgrounds presented above adequately 
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describe the data, and are in excellent agreement on the yields, polarization, and CP parameters of 
— > p^p^ decays. Statistical correlation between the two models is high (93% for the yields and /l, 
85% for and 87% for Cj^, as determined from a number of Monte Carlo experiments with event 
composition matched to the data). However, the two approaches differ in their sensitivity to the 
variations of the background composition, and have mostly uncorrelated systematic uncertainties 
associated with the B background PDFs and fit bias. Our two models represent the extreme 
approaches to describing the B backgrounds: one relies on the knowledge of the branching ratios 
associated with the dominant B backgrounds, and the other relies on the modeling of the kinematic 
distributions of B backgrounds in general. To further reduce these systematic effects, we average 
the results of the two models after the maximum likelihood fits. 

Other four-pion final states, such as B^ /o^vr+vr^ and B^ vr+Tr^Tr+vr^, require special care. 
The rates for these modes are not well constrained |22) . although their contributions are expected 
to be small in our invariant mass window. We parameterize the PDFs associated with these modes 
using the exclusive Monte Carlo samples, which assume uniform phase-space distributions of the 
final state mesons, and allow their yields to vary in the fit to the data. The absolute maximum of the 
likelihood occurs in the unphysical (negative) region for these yields, and we restrict the branching 
ratios for each of these modes to the range [0..6] x 10~^. The upper limit corresponds to several 
times the rate measured in Ref. |22) . With these restrictions, the yields for both B^ — > p^7r~^7r~ and 
B^ ir^n^TT^TT^ decays converge identically to zero in the maximum likelihood fit, indicating no 
evidence for any contribution from these modes to the fit region. 

5 RESULTS 

We find excellent agreement between our two approaches to modeling B background events. Ta- 
ble [1] shows the average results, while the difference between the two models is used to determine 
the systematic uncertainty associated with B background description. The B^ p^p^ decay is 
observed with a significance of 3.6a, as determined by the quantity \/—2 log(£o/^max)i where 
>Cmax is the maximum likelihood value, and Cq is the likelihood for a fit with the signal contri- 
bution set to zero. Both likelihoods include systematic uncertainties, which are assumed to be 
Gaussian-distributed. This significance level corresponds to a probability of background fluctu- 
ation to the observed signal yield of 2.8 x 10~^. We do not observe significant event yields for 
5° p°/o(980) and B° /o (980) /o (980) decays, nor of the non-resonant decays S° pOvr+vr- 
and B^ TT^-K^ir^ir^ . If the non-resonant contributions B^ p^tt^tt~ and B^ — > vr+Tr^vr+vr" 
were ignored in the fit, the significance for B^ — > p^p^ signal would go up to 4.2a. Background 
yields are found to be consistent with expectations. In Fig. [T] we show the projections of the fit 
results onto rriEs, Ai?, <S-shape, mi^2i and cos ^1^2 variables. Fig. [2] shows the distributions for the 
continuum qq component of the fit after subtracting the other components [23j . Fig. [3] shows the 
distribution of the likelihood ratio Lsig/ ^^Ci, where likelihoods Ci include all signal and back- 
ground PDFs. This ratio peaks near 1 for the signal events ( Fig. [3)3), and is highly peaked near 
zero for backgrounds. 

We also fit the proper-time distribution of the data sample, and determine the CP-violating 
parameters and for the longitudinal component of the B^ p^p^ sample. The results are 

listed in Tabled) The projection plots of At distributions for B^ and B^ tags, as well as the plot 
of the CP asymmetry, are shown in Fig. [H 
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Table 1: Summary of results: event yields (n), corrected for fit bias; fraction of longitudinal 
polarization (/l); selection efficiency (Eff) corresponding to measured polarization; branching frac- 
tion (Bsig), and significance including systematic uncertainties. The systematic errors are quoted 
last. We also show the background event yields for aiir, qq, charmless, and other BB components 
(statistical uncertainties only). 



Quantity Value 

n(SO ^ pV°) 85 ± 28 ± 17 

/l 0.70 ±0.14 ±0.05 

Eff (%) 23.8 ± 1.0 

Bsig (x 10-6) 0.84 ± 0.29 ± 0.17 

Significance, stat. only (a) 4.0 

Significance, syst. included (cr) 3.6 

0.4 ± 0.9 ± 0.2 

5f 0.5 ± 0.9 ± 0.2 

n{B^ /90/o(980)) -11 ± 16 

n(BO ^ /o(980)/o(980)) 6±6 

n{B° afiT^) 296 ± 42 

n( charmless) 348 ± 64 

n{BB) 2614 ± 134 

n{qq) 62298 ± 268 



6 SYSTEMATIC STUDIES 

The systematic uncertainties for all physics parameters are summarized in Table [2j Dominant 
systematic uncertainties in the fit originate from statistical errors in the PDF parameterizations, 
due to the limited number of events in the control samples, variations in the B background branching 
ratios fixed in the fit, and from the potential fit bias. The PDF parameters are varied by their 
respective uncertainties to derive the corresponding systematic errors. The fit bias is studied in 
a large number of Monte Carlo experiments, in which signal and charmless B background events 
are fully simulated, and b ^ c background events and continuum qq events are sampled from their 
respective PDF. The uncertainty associated with the B background model is 4 events for the signal 
yield, 0.01 for fi, 0.01 for C^^ and 0.11 for S^, derived from the difference between the fits for 
the two B background models. The systematic uncertainties due to the charmless background 
composition, arising from the uncertainties in the individual branching ratios and the CP content 
of B background [201 [2l] , are 5 events for the signal yield, 0.01 for Jl, 0.18 for and 0.14 for S'f . 
The above systematic uncertainties do not scale with event yield and are included in the calculation 
of the significance of the result. 
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We estimate the systematic uncertainty due to the interference between the -B — > p^p^ and 
a^TT^ decays using simulated samples in which the decay amplitudes for — > p'^ p^ are generated 
according to this measurement and those for B'^ afir^ correspond to a branching fraction of 
(33.2 lb 4.8) X 10~^ |25j . Their amplitudes are modeled with a Breit-Wigner function for all p — > vrvr 
and ai pn combinations and their relative phase is assumed to be constant across the phase 
space. The strong phases and CP content of the interfering state o]^7r^ are varied between zero 
and a maximum value using uniform prior distributions. We take the RMS variation of the average 
signal yield (14 events for the p^p^ yield and 0.03 for f^) as a systematic uncertainty. 

Uncertainties in the reconstruction efficiency arise from track finding and particle identification, 
and are determined by dedicated studies on copious control data control samples. Uncertainties 
due to other selection requirements, such as vertex probability, track multiplicity, and thrust angle, 
amount to 2.4% for the event yields, and are negligible for the polarization and CP observables. 



Table 2: Summary of systematic uncertainties. 



Source 


n{B 
fraction 


events 


fi 


5- 




Multiplicative 


Number of B mesons 


1.1% 










Event selection 


2.4% 










PID selection 


2.0% 










Track finding 


1.4% 










MC statistics 


<1% 




< 0.01 


< 0.01 


< 0.01 


aiTT interference 




14 


0.025 


0.07 


0.07 


Additive 


PDF variation 




6 


0.035 


0.07 


0.11 


Fit bias 




4 


0.014 


0.11 


0.09 


B background BR& CP 




5 


0.010 


0.14 


0.18 


B background model 




3 


0.005 


0.11 


0.02 


Total 


3.6% 


17 


0.047 


0.23 


0.24 



7 IMPLICATIONS FOR CKM ANGLE a 

To constrain the penguin contributions to B ^ pp decays, we perform an isospin analysis, by 
minimizing a term that includes the measured quantities expressed as the lengths of the sides 
of the isospin triangles. We use the measured branching fractions and fractions of longitudinal 
polarization of the i?^ P^P^ [6] decays, the measured branching fractions, polarization, and CP 
parameters and C^~ determined from the time evolution of the longitudinally polarized B^ — > 
p^p~ decay [7], and finally the branching fraction, the polarization, and the CP parameters and 
Cj^ of B^ p^p^ from this analysis. We assume uncertainties to be Gaussian and neglect 1=1 
isospin contributions, electroweak loop amplitudes, non-resonant, and isospin-breaking effects. 

Using the B^ p^p^ measurement we obtain a 68% (90%) CL limit on |a — aefr| < 14.5° 
(< 16.5°) where Oeff is constrained by the relation sin(2aefr) = '^'^"/(l ~ C*^^^)-*^/^. Fig. [5] shows 
the confidence level with or without using the measured CP parameters S^^ and C^^ in the isospin 
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analysis fit. We observe the four solutions around zero as in the B ^ tttt isospin analysis, but here, 
thanks to the additional constraint on S'j^, one of the four solutions a — = +11.3°. is favored. 

8 CONCLUSION 

In summary, we confirm our earlier evidence for — > p^p^ decays with a 3.6a" significance. We 
measure the p^p^ branching fraction of (0.84 it 0.29 it 0.17) x 10~^ and determine the 

longitudinal polarization fraction for these decays of /l = 0.70 it 0.14 it 0.05. We also constrain the 
CP parameters and S'£° for the longitudinal part of the S° —>■ p'^p^ final state: 

Sf = 0.5 ±0.9 ±0.2 
= 0.4 ± 0.9 ± 0.2 

where the first uncertainty is statistical and the second is systematic. These measurements com- 
bined with those for B'^ P^P^ and B^ — > p'^p~ decays provide a constraint on the penguin 
uncertainty in the determination of the CKM unitarity angle a. Wc find no significant evidence 
for the decays B° ^ pO/o, B° ^ /o/o, 5° ^ p^Tr+Tr" and 5° ^ Tr+ tt' Tr+ tt' . 
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Figure 1: Projections of the multidimensional fit onto (a) rriES, (b) AE, (c) event shape variable 
£, (d,e) di-pion invariant masses mi and 7712, and (f,g) cosines of the helicity angles cos0i^25 after a 
requirement on the signal-to-background probability ratio with the plotted variable excluded. This 
requirement enhances the fraction of signal events in the sample while keeping approximately 25% 
of signal events. The data points are overlaid by the solid blue line, which corresponds to the full 
PDF projection. The individual — > p^p^ PDF component is also shown with a dot-dashed red 
line. The sum of all other PDFs (including — > p^Jq and B^ /o/o components) is shown as 
the dashed black line. The D-meson veto causes the acceptance dip seen in (f,g). 
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Figure 2: Signal-subtracted projections (sPlots) of the multidimensional fit onto (a) m-ESj (b) AE, 
(c) event shape variable S, (d,e) di-pion invariant masses nii and m2, and (f,g) cosines of the helicity 
angles cos0i^2- The data are weighed to enhance the continuum qq events and effectively subtract 
all other fit components. The data points are overlaid by the solid blue line, which corresponds to 
the qq PDF component. 
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Figure 3: Likelihood ratio ^sig/ Yli^i> where hkehhoods Ci include all signal and background 
PDFs. Data points are overlaid hy a blue curve, which corresponds to the full PDF. The shaded 
teal histogram and the black dashed line correspond to the sum of background PDFs, and the red 
histogram corresponds to the signal contribution. Full range (a) and a zoom-in into the signal 
region (b) are shown. 




Figure 4: Projections of the multidimensional fit onto At variable for (top) B tags and (middle) 
tags. CP asymmetry Acp is shown in the bottom plot. The solid blue line represents the 
projection of the total PDF, the red dot-dashed line is the B^ p^p^ contribution, and the dashed 
black line corresponds to the sum of all backgrounds. A likelihood cut Cgig/ > 0.99 is applied 

to enhance the signal contribution. 
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Figure 5: Confidence level on a — aes obtained from the isospin analysis discussed in the text. The 
solid line CL includes the CP parameters Cj^ and in the fit. The dotted line corresponds to the 
usual isospin analysis without Sj^. The dashed curve is obtained without the two CP parameters 
Sj^ and Cj^. The horizontal dashed lines correspond to the 68% (top) and 90% (bottom) CL 
intervals. 
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